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1 Introduction 



There is a remarkable difference between the quantum properties of the clas- 
sically integrable and fully chaotic (ergodic) systems. While in the integrable 
case the wavefunctions possess an ordered structure, the cigcnstatcs of clas- 
sically ergodic systems appear random (Berry 1977, Voros 1979) in the semi- 
classical limit h ^ 0. They can be locally well represented as a superposition 
of plane waves with equal wavevector magnitude but random phases, leading 
to the Gaussian distribution of the wavefunction amplitude. The assumption 
of random phases may, however, break down on dynamical grounds when h is 
not sufficiently small. This leads to the phenomenon of scars, which are the 
regions of amplified wavefunction amplitude close to the short and weakly 
unstable classical periodic orbits (Heller 1984). 

In the classically mixed systems the energy (hyper) surface is split into 
both chaotic regions, within which the motion is ergodic, and the regular re- 
gions where motion is, as in the fully integrable case, quasiperiodic and con- 
fined to invariant tori within the energy surface. We call such systems generic 
since this is the most general type of Hamiltonian dynamics. The principle of 
uniform semiclassical condensation (PUSC, see Robnik 1988,1998) states that 
the quantal phase space distribution of any eigenstate, given by its Wigner 
function, should be uniformly distributed on a classically invariant object in 
the phase space when ^ — > 0. This object can be either an invariant torus in 
one of the regular regions in the phase space or a whole chaotic component. 

While the PUSC is valid for fully integrable and ergodic systems, its full 
potential is shown when applied to the generic case. Here it predicts that 
the eigenstates are separated into regular and irregular ones, depending on 
whether the classical object onto which their Wigner function condenses is 
a regular torus or a chaotic component, respectively. This has far reaching 
consequences leading, for example, to the picture of Berry and Robnik (1984) 
for the statistics of energy levels, where the regular and chaotic states are 
assumed to contribute independent level sequences to the total spectrum. 
This was confirmed by many numerical computations (Prosen and Robnik 
1993b, 1994 and 1999, Prosen 1995, 1996 and 1998, Robnik 1998). 

In this work we are interested in geometrical and statistical properties 
of high lying cigcnfunctions. It is an extension of previous work done on 
mixed type systems by Prosen and Robnik (1993c, 1994) and by Li and Rob- 
nik (1995a,b), but also relates to and further develops the paper by Li and 
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Robnik (1994) concerning the statistical properties of chaotic states, which 
drastically differ from the properties of regular states in classically integrable 
systems (Robnik and Veble 1998). The main step forward in the numerical 
direction was the use of the so-called scaling method, introduced by Vergini 
and Saraceno (1995), with which we obtained the states with consecutive 
indices around 2.5 ■ 10^ and thus it enabled us to go much farther into the 
semiclassical region. On the other hand, we were able to semiclassically re- 
construct the numerically obtained regular states by employing the Einstein- 
Brillouin-Keller (EBK) torus quantization in section 3.1, where we offer a 
new approach to this problem, especially on the numerical side. While we 
were not able to semiclassically reconstruct the chaotic states, we obtained a 
good prediction of their statistical properties by strictly employing the PUSC 
(see Robnik 1998) in section 3.2. 



2 Our catalogues of states 

We dealt with a model billiard system obtained by conformally mapping the 
unit circle with the complex quadratic polynomial, as introduced by Robnik 
(1983, 1984), 

z — > w{z) = z + , w{z^ = X + iy. (1) 

The range of x at y = inside the billiard is x G [—1 + A, +1 + A]. 
We used the value of A = 0.15, where the classical phase space is roughly 
equally divided into components of regular and chaotic motion. We chose 
the Poincare surface of section (SOS) to lie on the symmetry axis y = Q 
with coordinate x and the conjugate momentum px as the parameters of the 
surface. The intersection of the main chaotic component of our billiard with 
the SOS is shown in the figure 0. The coordinate x is taken relative to the 
center of the billiard (i.e. it is shifted by A w.r.t. a; = of equation (1), so 
that now the range of a; is a; G [—1, +1]), while is the x-component of the 
unit momentum vector. 

The quantum mechanics of billiards is described by the Helmholtz equa- 
tion 

(A + k^)i, = 0, (2) 
with the Dirichlet boundary conditions, where k'^ = 2mE/h'^. We limited 
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ourselves to the states with even parity with respect to reflection across the 
symmetry line y = 0. 

For each state we calculated the smoothed projection of the Wigner func- 
tion. The Wigner function of a state ip{c{) in the general case of degrees 
of freedom is defined in the full phase space (q, p) as (Berry 1983) 

W^(q,p) = ^^^ /c?^Xexp(-zp-X/;i)^t(q_x/2)^(q + X/2). (3) 

In our case the eigenfunctions il){x,y) generate their Wigner transforms 
W{x,y,Px,Py) through (|^), where N = 2. In order to compare the Wigner 
function of a state of our system with the classical SOS plot we took its value 
on the symmetry line {y = 0) and integrated it over py, 

psos{x,Px) = J dpyW{x,y = 0,Pa;,Py). (4) 

The result is 

Psos{x,p,) = ^J dXexp{-ip,X/h)i:\x-X/2,y = 0)ij{x + X/2,y = 0). 

(5) 

Here we see the reason for considering the even parity states only, because 
ip{x,y = 0) is exactly zero for odd states, and therefore a different approach 
must be used to analyze them. 

As is well known, the Wigner function is not positive definite but exhibits 
small oscillations that can blur the overall picture. We chose to smooth the 
projection of the Wigner function by a suitable Gaussian. It was chosen 
narrower than the minimum uncertainty Gaussian (in which case the Wigner 
function becomes the positive definite Husimi distribution) in order not to 
smooth out too many features, but still wide enough to reduce the oscilla- 
tions. 

The first catalogue of eigenstates and the corresponding smoothed Wigner 
function projections comprises the first 1000 even states. They were obtained 
by the conformal mapping diagonalization technique (as described in Robnik 
1984). We stress that by this method no levels and eigenstates were lost, as 
can happen with other approaches such as Heller's plane wave decomposition 
method (Heller 1984) and/or the boundary integral method (see e.g. Berry 
and Wilkinson 1984) etc. 
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Such a complete catalogue gives us a good picture of the overall behaviour 
of the system. Many of the states in this low energy region can be associated 
with the shortest classical periodic orbits. In a mixed type system such as 
ours, there are both stable orbits that are found within the islands of stability, 
or unstable ones that lie within one of the chaotic components. 

The shortest periodic orbits of our system are shown in figure These are 
the stable (labeled by 1) and unstable (5) periodic orbit with two bounces, 
the stable (2) and unstable (6) three bounce periodic orbit and the stable 
periodic orbit with four bounces (4). As our billiard is convex, there also 
exists an infinity of stable periodic orbits skipping along the boundary (3) 
of the billiard, which are associated with and support the whispering gallery 
modes (Lazutkin 1981,1991, Li and Robnik 1995a). 

Examples of the states that correspond to the stable periodic orbits with 
the indices (1-4) in the figure ^ are shown in the corresponding rows of the 
figure ^ with the smoothed projections of their Wigner functions shown in the 
figure |[ Each row shows the states of the same type with increasing energy. 
The states corresponding to the stable periodic orbits appear consistently and 
systematically across the catalogue of states. As we will see later, they can 
be attributed to the quantized tori in the regular regions, and can be found 
for those tori whose actions satisfy the Einstein-Brillouin-Keller quantization 
condition (see Robnik 1998). If we compare the Wigner plots of these states 
with the SOS plot in figure |l], we notice that the areas of greatest intensities 
of Wigner functions are found within the corresponding islands of stability 
in the SOS plot. 

The states that correspond to the unstable periodic orbits with indices 
(5) and (6) are shown in the rows of the figure ^ with the corresponding 
smoothed projections of their Wigner functions given in the figure with 
their energy again increasing along the rows. The states corresponding to the 
unstable periodic orbits emerge with varying intensities with respect to the 
background, and are less frequent with increasing energy. The approximate 
position of their emergence in the spectrum can be determined by the con- 
dition that the classical action along the periodic orbit should be a multiple 
of Planck's constant plus the Maslov's phase corrections due to the caustics 
formed by nearby trajectories (Robnik 1989). Therefore, in 2-dim billiards, 
the intervals (either in energy E or the cumulative number of states Af) of 
consecutive re-appearance of the eigenstates of the same type grow as the 
square root of E or Af. Such states with increased intensities of amplitude 
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close to the classically unstable periodic orbits are called scars (Heller 1984, 
1986, Bogomolny 1988). 

With increasing energy the longer periodic orbits start to manifest them- 
selves in the structure of eigenstates, however, many of the states are be- 
coming increasingly difficult to associate with simple periodic orbits. Two 
examples of such states are shown in figure 0. The left state is a regular state 
spanned by a torus in the neighbourhood of the stable periodic orbit with 
three bounces, characterised by two quantum numbers. It is interesting to 
note that while the quantum number along the direction of the periodic or- 
bit is quite large (about 140), the transversal quantum number is equal to 1, 
giving rise to a single nodal line along the direction of the periodic orbit. On 
the other hand the state on the right can not be associated only with a single 
unstable periodic orbit, but with a large portion of the chaotic component, 
this being more prominent in the smoothed projections of the corresponding 
Wigner function shown in the lower row of the figure |^ 

The separation of states into regular and irregular ones becomes fully 
explicit in our second catalogue of states. The catalogue consists of 100 
consecutive states starting at the consecutive index of about 2.5 ■ 10^. These 
states were obtained by the scaling method first introduced by Vergini and 
Saraceno (1995), that enables us to find a few states in the neighbourhood of 
a chosen wavenumber k. As this is a diagonalizational method no levels were 
missed. Almost each level in our small catalogue can be clearly identified as 
regular or irregular, an idea proposed already by Percival (1973). The only 
exception in the catalogue is a pair of states lying close together with respect 
to the mean level spacing, shown in figure §, where both of the states are 
superpositions of a regular {\ipr)) and irregular (lipi)) state. These two states 
are close to a degeneracy of two energetically equal but structurally different 
quantized classical objects. 

This mixing is the consequence of the fact that the regular and irregular 
state are not exact solutions of the Hamiltonian which leads to the matrix 
element 

H„={^r\H\iJ^) (6) 

not being equal to 0. When two such states, or, more precisely, their adiabat- 
ically corresponding states, are brought close together on the energy scale by 
varying a parameter of the system (e.g. A in our billiard), their eigenenergies 
do not cross but show a phenomenon of level repulsion (avoided crossing). 
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In the cases when the effect of other levels can be neglected, the smallest 
energy spacing between the two levels reached is twice the value of the ma- 
trix element (^). At this point the eigenstates are exactly the symmetric 
and antisymmetric superpositions of the well-separated states. If we vary 
the parameter of the system further, the adiabatic equivalents of the original 
states will exchange identities. 

While the avoided level crossings are typical for irregular states, the reg- 
ular states do not exhibit the phenomenon of level repulsion amongst them- 
selves (Berry 1983), except possibly on an exponentially small scale (due to 
the tunneling effects). Furthermore, according to the PUSC, there should be 
no level repulsion between the regular and irregular states as ^ — 0. As this 
repulsion is directly connected to the mixing of states, the relative number 
of mixed states, such as the pair shown in figure ||, is expected to tend to 
when effective Planck's constant of the system tends to 0. 

More states of both the regular and irregular type from this catalogue 
with their appropriate analysis will be shown in the next section. 

3 Analysis of states 

The main purpose of our work was to understand to what extent it is possible 
to describe the structure of individual eigenstates by semiclassical methods. 
For completely integrable systems there exists the Einstein-Brillouin-Keller 
method of torus quantization (see for example Berry 1983, Robnik 1998), 
where each semiclassical eigenstate is spanned by an invariant torus in clas- 
sical phase space, for which the classical actions are integer multiples of % 
with corrections due to the singularities of projection of the torus onto the 
configuration space. For fully chaotic systems, and general mixed systems 
(generic systems), the Gutzwiller periodic orbit theory (Gutzwiller 1990) can 
in principle be employed to obtain the density of states (the spectrum) and 
the wavefunctions, to the leading semiclassical approximation. 

The semiclassical methods cannot, however, predict individual energy 
levels within a vanishing fraction of the mean level spacing even in the limit 
^ — *^ (Prosen and Robnik 1993a, Robnik and Salasnich 1997). This limits 
their use in the analysis of statistical properties of spectra. Furthermore, in 
the chaotic case the levels are also very sensitive to perturbations of the sys- 
tem (Percival 1973). The Gutzwiller approach and method is very useful in 
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the qualitative analysis, and in certain context also quantitative (in describ- 
ing the collective and statistical properties), but it is still just the leading 
term in a certain semiclassical expansion, not good enough to resolve the fine 
structures with sufficient accuracy to make the analysis of individual states 
and energy levels reliable. Therefore, because of the sensitivity of chaotic 
cigcnstates and the approximating nature of the theory, to some extent the 
questions about the fine structure of individual chaotic states are irrelevant. 
It is therefore more appropriate to discuss the statistical properties of chaotic 
states, which are, however, less sensitive to perturbations. 

In a mixed system the phase space is divided into chaotic and regular 
components. Our work was guided by the principle of uniform semiclassical 
condensation (PUSC, see Robnik 1998), stating that when h tends to the 
Wigner function of any eigenstate uniformly condenses on an invariant object 
in phase space. This can be either a torus in the regular region or a whole 
chaotic component. Each state could thus be labeled as either regular or 
irregular (chaotic) in the semiclassical limit. By looking at the catalogue of 
states at high (and to some extent even at low) energies, one can see that this 
can indeed be done, though there is still the localization phenomenon present 
due to the still insufficiently low value of the effective Planck's constant. 

3.1 Regular states 

We start the analysis by considering the regular states. These are the states 
that can be attributed to quantized tori within the regular regions. For 
these states we tried to employ the EBK torus quantization. We construct a 
wavefunction on the torus as a sum of contributions 



of different projections j of the torus onto configuration space. Sj is the 
classical action with respect to some point on the torus and the classical 
density of trajectories on this projection. The phase of the wavefunction must 
change by an integer multiple of 27r when going around any closed contour 
of the torus. This gives us the quantization conditions 



(q) + 



(7) 
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where 7^ are the irreducible closed contours on the torus and rii the torus 
quantum numbers. The integers f3i are Maslov's corrections and arise due to 
the changes of phase 0j at the singularities of projection of the torus onto 
configuration space. At each caustic encountered along the contour 7j the 
wavefunction acquires a negative phase shift of 7r/2, and shifts by vr when 
reflected from a hard wallQ. From this consideration it follows that j3i counts 
the number of caustics plus twice the number of hard walls encountered along 
the contour. 

The main problem arises since, unlike in many completely integrable 
cases, the transformation to action-angle variables for regular components 
of mixed type systems is usually not known. The object we are dealing with 
is only the numerically calculated trajectory, so we must make the best of it. 

The task of finding the semiclassical EBK wavefunctions can be divided 
into two parts. The first one is finding the torus with the desired quan- 
tum numbers rii, the second one being the construction of its appropriate 
wavefunction in configuration space. 

In our billiard system all the regular tori except the whispering gallery 
ones wind around a stable periodic orbit with a finite number of bounces 
lb- We choose 61 to represent the movement of the trajectories along the 
corresponding periodic orbit while 62 represents the winding of the trajecto- 
ries around it. The Maslov index for the contour along 62 is = 2 since 
there are two singularities of projection (caustics) encountered, while for the 
contour along 9i this index is equal to Pi = 2lb. In the case of whispering 
gallery modes one similarly obtains /?i = and /?2 = 3, if 61 is taken along 
the boundary of the billiard and 62 'perpendicular' to it. 

We may find the appropriate torus by iteration. We start with a trajec- 
tory 7] in a regular island on the SOS and follow it sufficiently long until it 
returns to the desired neighbourhood of the initial point, thus approximately 
completing A^^i integer number of cycles in 61 and A^2 in ^2 in the time T. 
One can obtain the numbers Ni and N2 by knowing the total number of 
bounces A*";, and A^^^ of the caustics encountered during the process through 
the topological properties of the torus. For a torus winding around a periodic 
orbit, A^i = Nb/h and A^2 = Nc/2. The winding frequencies on the torus are 

^If the contour passes the singularity in the contrary direction to that of the Hamilto- 
nian flow on the torus, the phase shifts are of the opposite sign 
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then given by 



uji = 27rNi/T. 



(9) 



Finding the number of bounces Nb along the trajectory is straightforward. 
The number of caustics Nc may be obtained by starting a trajectory rj' near 
the original trajectory rj on the same torus (more are needed in case of 
more than two dimensions). Whenever the two trajectories cross there is 
a singularity of the density of trajectories and hence a caustic. 

One could of course use the monodromy matrix to find the caustics. There 
are two reasons for not doing so. The first one is the numerical simplicity of 
our approach. The more important reason is, however, that if the trajectories 
7] and f]' are taken too close together (infinitesimally close together in the 
monodromy matrix approach) we may observe caustics due to the possible 
graining of the desired torus into smaller islands of stability, which quantum 
mechanics is at the given value of effective Planck's constant still unable 
to resolve. A good criterion is that the trajectories should be separated by 
^ 1/nj for each 6i 

We still have to calculate the action integrals on the torus. The integral 
along 6*2 can be calculated by the integral 



along the curve 72 formed by the crossing points of the trajectory with the 
SOS. This integral is just l/27r of the area of the intersection of the torus 
with the SOS. Another action integral we can calculate is the action along 
the chosen orbit r], 



This integral is the sum of I = Nili + N2I2, so the integral along 62 is equal 
to 



We must iterate this procedure by choosing different starting points until 
the proper torus fulfilling the conditions has been found. In general there 
are as many conditions as there are degrees of freedom. For two dimensional 
billiard systems whose dynamics is independent of energy, all action integrals 
can be written in the form 




(10) 




(11) 



Il = {I-N2l2)/N,. 



(12) 




(13) 
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where Gm is an integral dependent purely upon the geometry of the classical 
object in question. In order to obtain the correct quantized torus, instead of 
solving two separate equations, due to this scaling property we need only to 
find the appropriate ratio of the (geometric) actions 

Ji ^ ni+/3i/4 

h n2 + /52/4" ^ ^ 

We used the robust bisection method to fulfill this condition since the depen- 
dence of actions upon initial conditions is not smooth due to the previously 
mentioned graining of the tori. Note that this procedure yields only the 
correct geometry of the quantized torus, the energy of which must still be 
determined by the quantization conditions 

The classical action as a function of time can then be written as 



/ L(q, q)dt = S{t) = {hui + hu:2)t, (15) 
J an 



rq(*) 
'qo 

so one immediately obtains the semiclassical energy as 



dS 

E = — = huJi + /2t^2- (16) 

In our case we took the quantized values for Ii and I2 and taken into account 
that for billiard systems the angular frequencies are of the form 



UJm = XmyE, (17) 

where Am are frequencies dependent again only upon the geometry of the 
chosen torus. 

In principle one could obtain the semiclassical energy without considering 
the angular frequencies Ui by simply finding the appropriate torus through 
the quantization conditions ^ and reading its energy. In a KAM system such 
as ours, however, these quantization conditions may be only approximately 
fuUfiled due to the fine structure of the phase space. In contrast to the actions 
Ii the frequencies Ui do not depend as strongly on the initial conditions (close 
to the main periodic orbit of an island of stability the transversal frequency 
typically does not vanish but is characteristic of the periodic orbit). So by 
taking the quantized values for actions Ii and the numerical values for uji in 
equation (|1^) the best estimate of the semiclassical energy is obtained. 
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Once the proper torus has been found, we proceed to the second step. 
We need to span the semiclassical wavefunction (J^) on this torus. Again 
we deal only with a trajectory, so we must find means of representing the 
wavefunction with it. We may write the wavefunction in the form 

^,(q) = ^lim 1/T £ dtSiqait) - q)Djiq,iit)) exp (^^Sjiq^iit) + 0,) , (18) 

where qci{t) is the classical trajectory in configuration space. This definition 
is appropriate only in the sense of integrals of the wavefunction over config- 
uration space. We can hence determine the function Dj by integrating the 
wavefunction over a small volume V, 

I, dVi^j = ^l™, (q) exp (^5, (q) + 0,) ^, (19) 

where T' is the time the trajectory spends in the volume V . This time is 
proportional to the density of trajectories A?j in the volume V , 

T' = TA^j{q)V'. (20) 

If we want the expression (|19D to be consistent with (|^), 

= (21) 

must hold. 

We still have to account for the phase shifts when the trajectory traverses 
from one projection of the torus onto another. This can be done hand in 
hand with the estimation of the density of trajectories Aj. Again we start 
a trajectory q^^ close to the original one (more trajectories are needed in 
more than two degrees of freedom). We could again use the monodromy 
matrix approach, but the same criticism applies here as in the case of finding 
the appropriate tori. Let us imagine a bundle of trajectories inside a small 
parallelogram spanned by the three points qd, q'^ and qd + qdSt, where all 
three points lie on the given torus. The area of the parallelogram is given by 
the absolute value of 

P = qd X (q'^, - qd)5t. (22) 
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As we are interested only in the relative sizes of this parallelogram, we may 
set 6t = 1. The reciprocal value of this area is proportional to the density of 
trajectories, 

A',=a/\P\. (23) 

The value of a does not change along the trajectory and can be subsequently 
determined by the normalization of the semiclassical wavefunction. When- 
ever the value of P changes its sign, the trajectory has encountered a caustic 
and has passed from one projection j of the torus to another one k. At this 
point the wavefunction acquires a phase shift 

0fc = 0i - vr/2. (24) 

The phase shift is equal to vr if a hard wall is encountered. 

The construction (0) is, in its original form, of course unsuitable for 
numerical computation. We must substitute the 6 function in the integral 
by a function of finite width. The other important reason for doing so is 
to again smooth the classical behaviour which can have a far more detailed 
structure than the quantum mechanics can yet resolve. This width may not 
necessarily need be isotropic. A good estimate for it is again to be of the 
order of l/rii in coordinates 6i projected onto configuration space. 

We used a Gaussian for the wide delta function. If we adjusted the width 
and the amplitude of the Gaussian so that it followed the classical density 
of trajectories, a remarkable similarity of our method with that of Heller's 
wavepacket approach was observed (Heller 1991). There are, however, two 
important differences. The wavepacket approach starts with a wavepacket 
with the expected value of energy equal to that of the exact wavefunction 
and then tries to construct its semiclassical approximation. Our approach is 
independent of the exact eigenstate. All the information we have to supply 
is the quantum numbers and the geometrical properties of the (projection 
of the) torus, obtaining from them both the semiclassical energy ( |I6D and 
the wavefunction. The other important difference is that the wavepacket ap- 
proach relies on the monodromy matrix of the trajectory, the use of which 
can be questionable to obtain semiclassical wavefunctions in mixed type sys- 
tems due to the fine structure of the classical phase space, as was already 
pointed out before. 

We show three examples of the regular states in the figures ^ to |Tl]. For 
each of the states we present the exact numerical quantum probability den- 
sity (top left), the probability density of its semiclassical approximation (top 
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right), the classical density of trajectories on the appropriate torus (bot- 
tom left) and the smoothed projection of the exact Wigner function (bottom 
right). The semiclassical wavefunctions shown are remarkable as they pos- 
sess all of the features of their exact counterparts that are larger than the 
appropriate wavelength. Note that for each torus there are two character- 
istic wavelengths since there are two quantum numbers associated with it. 
As it happens in our case, the two wavelengths can be of different orders of 
magnitude. 

There were, however, some 'regular' states that we were unable to semi- 
classically reproduce. The first class of these states can be described as local- 
ized chaotic states since their Wigner function clearly shows that they lie in 
the chaotic region, yet very close to an island of stability, which gives them a 
regular appearance. The other class are the states whose Wigner transforms 
lie in the regular regions, but where the primary tori have already been de- 
stroyed by the perturbation and now form secondary tori interwoven by small 
regions of chaotic motion, and even the smoothing of classical dynamics, as 
described above, fails. 

The accuracy of the semiclassical energies that we were indeed able to 
reproduce may seem remarkable, since the error is approximately 5 units 
of energy at the energies around 2 ■ 10''. Such accuracy, however, is still 
insufficient to perform short range spectral statistics since the mean level 
spacing in our system is approximately 8 units of energy. This experience is 
of course in agreement with the proposition and conclusion that the semiclas- 
sical methods (to the leading order) cannot resolve the energy spectra within 
the vanishing fraction of the mean level spacing, and also not the structures 
of the wavefunctions smaller than de Broglie wavelength (Prosen and Robnik 
1993a, Robnik and Salasnich 1997). 

3.2 Irregular states 

While for the regular states it was quite straightforward to find their semi- 
classical approximations, the nature of irregular states is very much different. 
The chaotic component of a system does not possess any obvious structure. 
While Gutzwiller's approach can yield the properties of a quantum system by 
a summation over all periodic orbits of its classical counterpart, the relevance 
of examining individual chaotic states becomes questionable. These states 
are very sensitive to small perturbations of the system, so in any physical sys- 
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tern the individual features of the states are lost when the effective Planck's 
constant tends to 0. The features that are insensitive to small perturbations 
are, however, the statistical properties of spectra and eigenstates. 

One measure of statistical properties of the wavefunctions is the wave- 
function autocorrelation function. 

The area of averaging e(q) close to the point q should be taken such that its 
linear size is many wavelengths across, however small enough that the local 
properties of classical mechanics within it are largely uniform. 

If one takes the Fourier transform of the Wigner function (^, it is easy 
to show that 

J W{q, p) exp(ip ■ x/h)d^p = V^^(q - x/2)V^(q + x/2). (26) 

By knowing the Wigner function of an eigenstate, it is then possible to use 
this result to calculate its autocorrelation function. 

According to the principle of uniform semiclassical condensation, the 
Wigner function of any chaotic state should uniformly condense on the whole 
chaotic component when the effective h tends to 0. Let us limit ourselves 
only to the cases of the Hamiltonians with an isotropic dependence upon p. 
We can write the semiclassical Wigner function in the form of a conditional 
delta function 

WnX^l, P) = «5({q, p} eV,-E- H{ci,p)) (27) 

where Pj denotes a chaotic component and a is the normalization constant. 
The Fourier transform (p6D of this Wigner function is 

Wt),{^, p)p^"^ ey:p{ip-x/h)dpdnp = a-^ — - / dVlpexpiip-x/h), 

If (q,p(q)) J^v^-^M) 

(28) 

where p(q) denotes the absolute value of momentum at the point q. The 
integration over the spatial angle VLp is performed along all the directions of 
momentum that constitute the chaotic component T>i at the point q. The 
autocorrelation function is then equal to 

^ I . _ (/Q,G2?.(q')^^peXp(zp-x/;i))q,e,(q) 

ViflpeACq') ""p/q'ee{q) 
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The averaging area should again stretch across many wavelengths. 

If the chaotic component is equal to the whole energy surface, as is the 
case in completely ergodic systems, in the case of two degrees of freedom one 
obtains the well known Berry's result (Berry 1977) 

Cergodic(x) = Mpici)r/h), T = |x|. (30) 

However, when the system is of the mixed type the autocorrelation function 
ceases to be isotropic and acquires contributions of higher order Bessel func- 
tions. We can obtain these contributions by rewriting integrals /^^gp. f{4>p)d(f)p 
by integrals of the characteristic function / XVi{4>p) f {<f>p)d4>p, where in two 
degrees of freedom the spatial angle is replaced by a simple angle (pp and 
f{4>p) is an arbitrary function of (f)p. If we write the characteristic function 
as a Fourier series, 

oo 

Xi?i(q;0p)= J2 '«m(q)exp(im0p), (31) 

m=— oo 

it is quite straightforward to show by using the integral representations of 
the Bessel functions that 



(Em=-oo ^m(qO^ "'^m(p(qOr/;i)e xp(^m0^))q/ee(q) 

(«o''(q'))q'Gc(q) 



r< X _ \Z^m=-oo '^m o-mV/^V^;' /"';"=^PV^"^rx;/q'Ge(q) /go\ 



where c^^- is the polar angle of the vector x. 

As in the case of the regular states, numerically we cannot deal with com- 
ponents of phase space but with trajectories. So we start a trajectory within 
the chaotic component Di and not the direction of its momentum 0^(q) at 
each passage j through the neighbourhood of the point q. The averaged 
characteristic function for this neighbourhood can then be represented as 

xx)iiq, 0pj/'«o iqj ^ — > v>^) 

where dj are the lengths covered by the particle in the averaging neighbour- 
hood at the passage j, and n is the number of passages. One can check that 
this method tends to the proper angular distribution as n — > cxo for any 
shape of the averaging neighbourhood, if one assumes the homogeneity of tra- 
jectories at a given angle (which is true for a small enough neighbourhood. 
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over which the classical phase space picture does not vary). For any incident 
angle the conditional expected value of a contribution to the equation (^) 
is proportional to / djdx, where x is the homogeneously distributed impact 
parameter (the direction perpendicular to the incident angle). This integral 
gives just the area of the averaging neighbourhood and is clearly the same 
for all incident angles. 

If the area of averaging is large enough so that the variations in the 
classical phase space picture become important, the above procedure is still 
valid as long as the value of the momentum does not change appreciably (as 
is the case in our billiard system, where between the bounces the momentum 
remains constant). One can then imagine the large averaging area as being 
cut into smaller ones within which the above assumptions still hold true. 

We compared the autocorrelation functions for a few chaotic states with 
the semiclassical prediction in figures O to The averaging area e(q) 



was taken as a circle of radius 0.2 around the point {x,y) = (0.65,0) (the 
coordinates are as defined in equation (1)). It was taken the same for both 
the semiclassical prediction and for the numerical results. The averaging 
radius was taken quite large in order to reduce the localization properties of 
the wavef unctions, which are still apparent at the values of effective h that 
we were able to obtain. But this radius still has to be taken small enough in 
order not to completely smooth out the classical dynamics. The agreement 
with the semiclassical prediction is quite good particularly in figure 14. In all 
cases it clearly deviates from the Berry's prediction for fully ergodic systems 
(|30D, as it must for mixed systems, and tends towards our semiclassical result. 
Although in some cases there are amplitude deviations from our prediction, 
in all of the plots the phase of the numerical correlation function matches 
the phase of its semiclassical prediction and is significantly different from the 
phase predicted for the fully ergodic case. 



4 Discussion and conclusion 

As already presented by Prosen and Robnik (1993c), the classification of 
states into regular and irregular ones is well founded when the effective 
Planck's constant tends to 0. Its theoretical foundation is the principle of 
uniform semiclassical condensation of Wigner functions of eigenstates (Rob- 
nik 1988, 1998). This separation is not strictly a semiclassical phenomenon 
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since even some of the lowest levels in our catalogue of low lying states can be 
classified as either regular or irregular. In the high energy catalogue of states 
each state can easily be classified as either chaotic or regular, with only one 
notable exception, where two close lying states are a superposition of a reg- 
ular and an irregular state. These exceptions are expected to disappear with 
higher energies. While the states can be separated with respect to classical 
dynamics, the chaotic states in our high energy catalogue still exhibit the 
phenomenon of dynamical quantal localization. Their Wigner functions are 
not uniformly extended over the whole chaotic component, but are significant 
only on a part of it. This localization is expected to disappear at sufficiently 
small effective fi, when the quantum mechanical break time tbreafc = h/ AE, 
where AE is the mean level spacing, becomes longer than the time for a 
typical trajectory to explore the whole chaotic component (diffusion time). 

The most important part of this work is of course the semiclassical analy- 
sis of states. We were able to reconstruct both the semiclassical wavefunction 
and the semiclassical energy of the regular states by using the EBK quan- 
tization. Since our system is of a KAM type, most of the resonant tori are 
destroyed, forming smaller islands of stability interwoven with chaotic com- 
ponents. The classical mechanics thus shows a rich structure that quantum 
mechanics at a fixed value of h is still unable to resolve. In order to obtain 
the regular semiclassical wavef unctions, we had to appropriately smooth out 
this fine classical behaviour as explained in the section |3.1| . 

We were unable to predict the individual properties of chaotic states. We 
made a step forward, however, in describing their semiclassical statistical 
properties. We obtained a semiclassical prediction for the autocorrelation 
function of their wavefunction, which differs from the one for fully ergodic 
systems as it is not isotropic. The numerical results confirm this prediction, 
although there are still localization phenomena at the currently attainable 
effective % that cause deviations from it. 

One aspect that needs to be investigated further is the localization prop- 
erties of the chaotic states. As Casati and Prosen (1998) show in the example 
of a fully chaotic stadium billiard, in the diffusive regime (e-stadium, having 
very large ergodic time), the quantum diffusion is stopped by the cantori in 
phase space leading to localization. How these and similar ideas translate to 
the case of a mixed type system remains so far an open question. 
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Figure 1: The SOS section of the main chaotic component of the A = 0.15 
bilhard. The coordinate x here is shifted by A to the right so that x G 
[— 1, +1], and X = corresponds to x = A of equation (1). 

Figure 2: The shortest periodic orbits of the A = 0.15 bilhard (explained 
in text). The stable orbits are shown with full and the unstable ones with 
dashed lines. 

Figure 3: A selection of eigenstates corresponding to the stable periodic 
orbits 1, 2, 3, 4 in the figure Q from top to bottom, respectively. In the top 
half of each plot we show the eight equally spaced contours of the probability 
density from to its maximum value. In the bottom half we show the nodal 
lines of each wavefunction. 

Figure 4: A selection of eigenstates corresponding to the stable periodic 
orbits 5 and 6 in the figure H in the top and bottom row, respectively. In 
the top half of each plot we show the eight equally spaced contours of the 
probability density from to its maximum value. In the bottom half we show 
the nodal lines of each wavefunction. 

Figure 5: The corresponding smoothed projections of the Wigner functions 
for the wavefunctions in the figure |[ The contours are spaced in ten intervals 
from to the maximal value. The negative value contours are plotted at the 
same spacing but with thinner lines. We do not plot the zero level contour 
due to many oscillations of the Wigner functions when it is close to zero. As 
in all subsequent Wigner plots, the momentum is measured in units of p', 
which is related to the energy of the eigenstate hj E = p'"^. The locations 
of the relevant classical periodic orbits are marked by the "bullets". The 
coordinate x here is shifted by A as in figure 1. 

Figure 6: The corresponding smoothed projections of the Wigner functions 
for the wavefunctions in the figure ^. The plotting method is the same as in 
figure ^. 
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Figure 7: An example of the probability density for a regular (top left) and 
an irregular (top right) state with the corresponding smoothed projections of 
their Wigner functions shown below, with the same plotting methods as used 
in corresponding previous plots. In the phase space plots (bottom row) the 
coordinate x is shifted by A as in figure 1. The "bullets" mark the location 
of the stable period 3 periodic orbit, which is the skeleton of the quantized 
invariant torus. 



Figure 8: The probability density (16 equally spaced contours from to the 
maximum value) for a pair of close lying states (fc^ = 20421106.7347 left 
and k"^ = 20421107.0691 right) that are superpositions of a regular and an 
irregular state. 



Figure 9: The probability density for a regular state with k"^ = 20420831.0603 
(top left), its semiclassical approximation with fc^^ = 20420828.18 (top right), 
with the quantum numbers on the torus being rii = 6752 and n2 = 37 (16 
equally spaced contours from to the maximum value in both cases). In 
the bottom row we show the classical density with 20 contours from to the 
maximum of the appropriate torus (left) and the smoothed projection of the 
exact Wigner function (right), with ten contours from to the maximum 
value. In the phase space plot (bottom right) the coordinate x is shifted by 
A as in figure 1. 



Figure 10: The same as in figure | but for the state with P = 20421002.7443 
and its semiclassical approximation with k"^^ = 20420999.46, rii = 2861 and 
ns = 20. 



Figure 11: The same as in figure |^ but for the state with k"^ = 20421387.1741 
and its semiclassical approximation with k^^ = 20421385.96, rii = 8648 and 
n2 = 1. 
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Figure 12: In the top row we show the probabihty density for the chaotic state 
with k"^ = 20420756.1273 (left, 8 contours) with the smoothed projection of 
its Wigner function (right, 10 contours). The circle of radius 0.2 centered 
at X = 0.5 is the region of averaging. In the phase space plot (top right) 
the coordinate x is shifted by A as in figure 1. In the bottom row we plot 
the wavefunction autocorrelation function (averaged over the small circle as 
explained in text) in the x (left) and y (right) directions. 



Figure 13: The same as in figure but for the state with k"^ 
20421005.3834. 



Figure 14: The same as in figure but for the state with k'^ 
20421262.6667. 
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